use "$Data_generated/Commune.dta", clear
gen Pasc_bin=0 if Pasc<0.3365316 & !missing(Pasc)
replace Pasc_bin=1 if Pasc>=0.3365316 & !missing(Pasc)
replace Pasc_bin=. if missing(Pasc)
cibar nb_sectperm if year==1926, over(Pasc_bin) barc(gs12 gs3) graphopts(ytitle("Avg Nb of Branches") legend(order(1 "Relig<Median" 2 "Relig>Median")))
graph save "graph1.gph", replace
cibar nb_sectperm if year==1935, over(Pasc_bin) barc(gs12 gs3) graphopts(ytitle("Avg Nb of Branches") legend(order(1 "Relig<Median" 2 "Relig>Median")))
graph save "graph2.gph", replace
graph combine "graph1.gph"  "graph2.gph", ycommon
graph export "$Figures/Figure4b.pdf", name(Graph) replace
graph export "$Figures/Figure4b.png", name(Graph) replace